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■ We study the long-time asymptotic of the solutions to Maxwell’s equation in the case of a upper- hybrid 
resonance in the cold plasma model. We base our analysis in the transfer to the time domain of the recent 


results of B. Despres, L.M. Imbert-Gerard and R. Weder, J. Math. Pures Appl. 101 ( 2014) 623-659, where 
,_l the singular solutions to Maxwell's equations in the frequency domain were constructed by means of a 
limiting absorption principle and a formula for the heating of the plasma in the limit of vanishing collision 
frequency was obtained. Currently there is considerable interest in these problems, in particular, because 

S upper ' hybrld resonance8 are a po88lble 8cenart ° for the heatlng of pla8ma8 ' and 8ln “ they can be a ™ del 

tO for the diagnostics involving wave scattering in plasmas. 
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1 Introduction 


In plasma physics (5j (7, 17) upper-hybrid resonances may develop at places where a density gradient of 
charged particles excited by a strong background magnetic field generates singular solutions to Maxwell’s 

equations. This phenomenon shows up in propagation of electromagnetic waves in the outer region of the 
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atmosphere, as explained first in [6]. It also appears in reflectometry experiments 1101 [Til and in heating 
devices in fusion plasma |16l in Tokamaks. An important feature in this direction is the energy deposit 
which is finite. It may exceed by far the energy exchange which occurs in Landau damping 1171 1251 . 
Notice however that there exists situations where hybrid resonance and Landau damping are modeled in 
a unique set 11311551 . Furthermore, our model could be applied in diagnostics involving wave scattering 
at the upper-hybrid resonance l2ll . 1 101 . 


The starting point of the analysis is the linearized Vlasov-Maxwell's equations of a non-homogeneous 
plasma around a bulk magnetic field B 0 ^ 0. It yields the non-stationary Maxwell’s equations with a 
linear current 

{ -^iE + VAB = /joJ, J = -eJV e u e , 

( 1 . 1 ) 


~~ jr<9tE + V A B — fx 0 J, 

<9 t B + V A E = 0, 

m e d t u e = -e (E + u e A B 0 ) - m e v\i e . 


The electric field is E and the magnetic field is B. The modulus of the background magnetic field |B 0 | and 
its direction b 0 = will be assumed constant in space for simplicity in our work. The total magnetic 
field is expanded as first order as B to t = B 0 + B. Note that in the last equation in 11.11 B is neglected. 
The absolute value of the charge of electrons is e, the mass of electrons is ro e , the velocity of fight is 


c = yj where the permittivity of vacuum is e 0 and the permeability of vacuum is /i 0 . The third equation 
corresponds to moving electrons with velocity u e , and the electronic density N e is a given function of the 
space variables. One assumes the existence of a bath of particles which is the reason of the friction 
between the electrons and the bath of particles with collision frequency u. Much more material about 
such models can be found in classical physical textbooks (5J [171. The loss of energy in domain <> can 
easily be computed in the time domain starting from 11.11 . One obtains 

12 


d 

dt 


eo|E[ 


|B| 2 . m e N e \u e \ 2 


2p 0 


= — / vm e N e |u,, | + boundary terms. 


Therefore, 


Q(v) = / vm e N e |u e 
J n 


( 1 . 2 ) 


represents the total loss of energy of the electromagnetic field plus the electrons in function of the collision 
frequency v. Since the energy loss is necessarily equal to what is gained by the bath of particles, it will 
be referred to as the heating. In fusion plasma, a value of v s=s 10 -7 in relative units is often encountered. 
It is therefore tempting to set the friction parameter, i.e. the collision frequency v , equal to zero, but this 
naive approach is incorrect, as we explain below. 

Equations 11.11 can be written in the frequency domain, where w is the frequency, that is d t = —ico 
where for simplicity of the notations ui > 0. We assume that the bulk magnetic field is along the 2 
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coordinate, i.e. b 0 = (0,0,1). We obtain. 


E + V A B = —fj, 0 eN e u e , 

—*wB + V A E = 0, 

—im e um e = — e (E + u e A Bq) — m e v\i e . 


(1.3) 


One can compute the velocity using the third equation wu e + u c iu e A b 0 = —— iE where the cyclotron 


frequency is w c = 


e|B 0 


■. The frequency u = w + iv is shifted in the complex plane by a factor equal to the 


friction parameter. It is then easy to eliminate u e from the first equation of the system 11.31 and to obtain 
the time harmonic Maxwell’s equation 
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V A V A E - 


0s< 


idE = 0. 


(1.4) 


The dielectric tensor is the one of the cold plasma approximation, the so-called Stix tensor, [7] fl7 1 
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(1.5) 


The parameters of the dielectric tensor are the cyclotron frequency w c and the plasma frequency u p = 
J which depends on the electronic density N e . We are interested in the physical situation where the 
electronic density N e is not constant, that is V N,, ^ 0. Observe that the cyclotron frequency w = w c is 
a singularity of the dielectric tensor. In this paper we consider u ± uj c , hence, the dielectric tensor 11.51 
is smooth. In fact, we focus on the paradoxical upper-hybrid resonance that appears when w = iv h := 


jjl + In plasma physics coh is called the upper-hybrid frequency. 

2 _ 2 

If we set v = 0 the first two diagonal entries in the Stix tensor are equal to . The crucial issue is 

that they are equal to zero when u = ujh and that they continuously change in sign when w increases from 
values smaller than uj h to values bigger than uju . For this reason (see [15)) the system 11.41 with u = 0 is ill 
posed: the solution is not unique and, furthermore, it has singular solutions that contain distributions. 
The way out of this dilemma ll5l is a limiting absorption principle, we take v to zero in a limiting sense: 
We consider v > 0 small and we construct an unique solution to 11.41 characterized by its behaviour at 
the point in space where ujh = w and by demanding that it goes to zero at spacial infinity, away from 
the sources of the electromagnetic field (the system 11.31 is assumed to be coercive (non propagating) at 
infinity). We call it the singular solution. We then prove that as v l 0 the singular solution converges (in 
distribution sense) to a limiting singular solution that is the appropriate physical solution to our problem. 
Furthermore, we give a formula for the limiting heating Q(0 + ) := lim^o Q(V), that turns out to be positive. 
The fact Q(0 + ) > 0 implies that the hybrid resonance is able to transfer a finite amount of energy from 
the electromagnetic field and the electrons to the bath of particles (i.e. to heat the bath of particles) 
even in the limit when v = 0. This is, indeed, a remarkable result. The physical interpretation is that 
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as v goes to zero -and the solution becomes singular- the velocity of the electrons increases (since the 
friction with the ions goes to zero) and there is a compensation in the right-hand side of 11.2) . In the end, 
the increase in the velocity of the electrons dominates, so that, in the limit the heating <2(0 + ) is positive. 
Moreover, from the mathematical point of view, the fact that the singular solution with v > 0 converges to 
a limiting singular solution implies that for v small the singular solution is close to the limiting singular 
solution and, hence, it does not change very much with v. This means that a small v positive can be 
used as a regularization parameter to numerically compute the limiting singular solution and the limiting 
heating Q(0+). This numerical scheme has been successfully used in the Ph. D. thesis |20| that contains 
extensive numerical calculations. It was found that the numerical solution with small v > 0 converges 
fast to the exact solution (with v = 0) in point-wise sense, except of course, at the singularity. Moreover, 
a large fraction of the energy of the incoming wave may be absorbed by the bath of particles, up to 95% in 
the case of normal incidence, and up to 76.7% in the case of oblique incidence. Our results in jl51 , |20| , in 
particular our formula for the heating Q(0 + ) shows, in a rigorous and quantitative way, that upper-hybrid 
resonances are, indeed, an efficient method to heat the bath of particles. 


We now present, for later use, our results in |15| in a precise way. To study the upper-hybrid resonance 
we consider the 2x2 upper-left block in 11.41 . that corresponds to the transverse electric (TE) mode, 
E = (E x ,E y , 0), where the electric field is transverse to the bulk magnetic field B 0 . We assume that we 
have a slab geometry: all coefficients in )1.5) depend only on the coordinate x. Furthermore, we suppose 
that that E x ,E y , are independent of 2 , that is the coordinate along the bulk magnetic field B 0 . 


Then, in the limit case v = 0, the 2x2 upper-left block in 11.4) gives 


w 

~\~d y E x 

-d x Ey 

= 0 , 


dyW 


—iSE y 

= 0 , 

( 1 . 6 ) 

d x w 

-\- iSE x 

—aE y 

= 0 , 



where we find convenient to introduce the vorticity W := d x E y — d y E x that is proportional to the magnetic 
field B z . The coefficients a, 6 are equal to 




and 


0 = — X 




w (w 2 - w 2 )' 


(1.7) 


In plasma physics the system [L6l is called the equations for the X-mode, or the extraordinary mode, and 
also the extraordinary wave. 


We consider the system 11.61 in the two dimensional domain, 


f l = {(x,y) € R 2 , — L < x, y £ R, L > 0} . 


We assume the non-homogeneous boundary condition 


W + i\n x E y = g on the left boundary x = —L, A > 0, 


( 1 . 8 ) 
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that corresponds to a source, like a radiation antenna that is used to heat the plasma.. We suppose that 
a and S satisfy conditions that correspond to a upper-hybrid resonance at x = 0. The main assumptions 
are: The function a is twice continuous differentiable and 6 is continuous with continuous first derivative. 
Furthermore, a(0) = 0, a'(0) < 0, and a ^ 0 for x / 0. Furthermore, 5(0) ^ 0. We assume that the 
coefficients are constant far enough from the source, i.e. : that there exist 5^, , and H > 0 such that 

S(x) = 5oo and a(x) = H < x < oo, 

and that our system is coercive (non propagative) at infinity, 

<£,-£,>0. (1.9) 

In 1151 we also assume other technical conditions. For a complete list see assumptions Hi — H 6 in 1151 . 

For the purpose of considering the limiting absorption principle we actually consider a model where the 
Stix tensor 11.51 is replaced by e(v) = e(0) + ivl , where, of course, e(0) is the value of the Stix tensor 11.51 
for v = 0. This is a simplified linear approximation. For the two by two block in 11.51 . that corresponds 
to the upper-hybrid resonance, this amounts to replace a by a + iv in 11.61 . We denote the solution 
by, (/I)), Ky, IT !/ ). Note that when v > 0 our system 11.71 is well posed. Clearly, v in our linear model 
is not exactly the same quantify as the collision frequency v in 11,11 11.21 11.31 11.41 . However, since it 
plays the same role, we still denote it by v and we call it the collision frequency, to simplify the notation. 
Furthermore, we take advantage of the fact that the coefficients a , S only depend on x to take a Fourier 
transform along y, 

( \ ( El \ 

u e ’ ! ' = V e ’ v := / e~ i6y Ey dy. 

\ W 6 ’ 1 ' ) , ' R \ W v ) 

Then, \J llM is a solution to the system 

( W 0 ' v + i6U e ’ v - ^V e ' v = 0, 

J iOW 0 ’’' -(a{x) + iv)u e ' u -i5{x)V e ^ = 0, (1.10) 

l -£ w8 ' U + tiix)U 0 ’ v - (a(x) + iv)V e ' v = 0. 

We denote by 

\l 0 ’ v = (U 0 ^,v 0 ’ v ,w 9 ^) 

the solution to II. 101 that is uniquely defined by the following two conditions [ 15]: 

1. It is exponentially decreasing at infinity. Recall that by 11.91 the system II. 101 is coercive (non 
propagative) for x > H. 

2. Its value at the origin is normalized with the requirement 

ivU 6 ’ v { 0) = 1. (1.11) 
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Note that as v 0 the value of U°■''{{)) goes to infinity and for this reason we call it the singular solution. 
Observe that in 051 we denote this solution by U®’ 1 ' to distinguish it from other solutions introduced 
there. 

In Proposition 5.16 of 1151 we proved that. 


limU e,I/ = U e ’ + := f P.V. , 
t>lo \ a(x) cP(0) 

where u e,+ ,v e,+ ,w 0,+ € L 2 (—L, oo) and 8d is the Dirac delta function at the origin. 


1 


H- 777 - 7 S D +U e - + , v e ’ + , W 2 + 


( 1 . 12 ) 


The limit U e + is a solution in the sense of distributions of the system II. 101 with v = 0. We call it the 
limiting singular solution. Furthermore, in Remark 10 of 1151 we prove that the functions u 2 + ,v 2 + ,w e 2 + 
can be obtained as the unique, properly normalized, solution of the system of equations. 


w. 


e,+ 


,e,+ 


idu, 


e,+ 


i6w 2 + — a(x)u 2 ’ + — iS( x)i 


h+ 


= -id P.V. -f, 

a[x) 


~-t w 2 + + i$(x)u 2 ’ + - a(x) 


e,+ 


- 


= -iP.V. 


s ( x ) 
a(x ) 


On r 
a'( 0 ) 

5 ( 0 ) 71 - c 
7m 6 D- 


There is a striking similarity with scattering theory 1311 . where the solutions obtained by the limiting 
absorption principle are characterized as the unique solutions that satisfy the radiation condition, that is 
to say, they are uniquely determined by the behavior at infinity. Here, in order that the singular solutions 


are uniquely determined we have to specify their behavior at +oo and also their singular part at x = 0. 
P. V. + jpjf-rdi). From the mathematical point of view it is natural that we have to specify the singular 
part because our equations are degenerate at x = 0. The physical interpretation of this fact is that x = 0 


acts as a scattering channel, and then, a condition has to be specified at x = 0. Similar situations arise 
in scattering theory, in the case of obstacles with a singular boundary 1181 and for manifolds with ends 

m. 


We construct the physical singular solution (in distribution sense) to 11.61 by inverting the Fourier 
transform along y and by demanding that the boundary condition 11.81 be satisfied. Let us designate by 
g the Fourier transform of the function g that appears in the boundary condition 11.81 . 

W) ■= [ g{y)e~ l6v dy. 

JR 


In Theorem 1.1 of fra we prove the following results: If g e L 2 (R) with g of compact support, there 
exists a solution (in distribution sense) of 11.61 with boundary condition 11.81 that goes to zero infinity. 
This solution, that we call the physical singular solution, has the following representation formula 

( X , y ) = -L [ ^lv e ^ e if) yde. (i.i3) 

2tT Jr T e - + 
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The transfer coefficient r 6 ' + is given by, 


T 0,+ := w e ’ + {-L) + i\V e ’+{-L). 


The electric field E x does not belong to oc ((—L,oo) x R) (of course, unless the source g is identically 
zero). The other components are always more regular: In particular E+,W + e L 2 ((—L, oo) x M), but they 
have a logarithmic singularity at x = 0, see equations (5.57), (5.58) in [151 . 

Moreover, the value of the resonant heating is 



(1.14) 


Note that it is known since the seminal work of Budden go, m that the solution to the X-mode 


equations (in a particular case that has explicit solution ) have a singularity as 1 /x at the position (x=0) 


of the upper-hybrid resonance. In the recent paper Il4l the singularity around the hybrid resonance has 
been studied with a purely local method based on deformation in the complex plane. Our results and 
our method in 1 1 51 are fundamentally different. Note that the fact that a singularity as l/x ( that is not 
integrable in a neighborhood of zero ) appears in the solution means that it is not a solution in distribution 
sense. Of course, the term l/x can be interpreted as a principal value, but also a Dirac delta function, 
or another distribution with support at x = 0 could be added. Our global analysis in 1151 , based in the 
limiting absorption principle and the requirement that the solution goes to zero at infinity, is the crucial 
fact that allowed us to calculate the singularity of the physical solution given in 11.13) . Furthermore, this 
was essential to compute the heating )1.141 that depends on the singularity of the physical solution. To 
our knowledge our paper f 1 51 is the first that gives a formula for the heating of an upper-hybrid resonance 
with a discussion to the radiation condition at infinity. 

The results of IT51 pose challenging new mathematical and physical questions, for example: 

OPEN PROBLEM 1.1. To study the propagation, in the time domain, of a wave packet (a definition in the 
context of waves in plasmas is to be found in (5], see also l23l 1271 ) that is incident from the left on the 
upper-hybrid resonance. 

Our results in 11 51 in the frequency domain make it intuitively clear to expect that part of the incoming 
wave packet will be reflected and that the transmitted part will be absorbed by the bath of particles as 
heating. To set-up this translation of our results from the frequency domain to the time-domain is an 
issue that has to be treated with mathematical care because of the singularity of the physical solution. In 
Section 2 we solve this problem. We prove that as i -> ±oo the part of the wave packet in [— R, —L] tends 
to zero for every R > 0. Since there is no heating for x < — L this implies that the reflected part of the wave 
packet actually travels to —oo in the x direction. Furthermore, we prove that as t —> ±oo the wave packet 
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tends to zero in the domain x > L. Since we assume that our problem is coercive (non propagative) for 
x > H the wave packet can not travel to plus infinite in the x direction. This implies that for large times 
the transmitted part of the wave packet is absorbed by the bath of particles as heating. 


OPEN PROBLEM 1.2. Determine a general perturbation theory and appropriate limiting absorption and 
limit amplitude principles for the Maxwell equations with hermitian complex valued continuous dielectric 
tensors e(x) = e(x)* with continuous eigenvalues which change sign. 

Here perturbation is for the dielectric tensor. For example one may use e = e-\-vI 3 with I 3 the identity 
matrix in C 3 , or more generally, as in the case of the the cold plasma dielectric tensor 11.51 which admits 
the expansion = e Q + ?D + 0(v 2 ) with D a physically based hermitian non negative matrix that depends 
linearly on v 1141 . Proving a (generalized) limiting absorption principle in these situations is an important 
and challenging issue. In this direction some inspiration can perhaps be obtained by comparison with 
other non-standard dielectric tensors and related physical and mathematical problems that can be found 
in |2l[3U4U8l|29|l30l. 

Roughly speaking, the limit amplitude principle means that in the case of Maxwell equations in the 
time domain with a periodic in time force term, the solutions oscillate with the frequency of the force 
for large times. For standard dielectric tensors the limit amplitude principle follows from the limiting 
absorption principle Il2[ 1221, 28j . It is important to extend the validity of the limit amplitude principle to 
non-standard dielectric tensors as above. 

These are, indeed, challenging problems whose solution does not follow from the standard results in 


fll 151H21 1221 1251 1251. 

2 Large time asymptotic of the solution to Maxwell’s equations 


To study the propagation of a wave packet we consider the problem in the whole space. We suppose that 
the assumptions that we made in the Introduction and in 1151 are valid for x > —L, and furthermore, that 
N e (x) = 0 for x < —L. Then, for —oo < x < —L we have an isotropic and homogeneous medium (vacuum) 
with permittivity e 0 and susceptibility // 0 . Hence, for x < L and after taking a Fourier transform in the y 
variable we see that U 0 + must be a solution to the system. 



( 2 . 1 ) 


for x < —L. 
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If we eliminate U e - + in 12.11 we obtain the following system of ordinary differential equations 


where, 


d ( V e ’+ 
di ^ W e '+ 


A e {x) 


y8,+ 

W 9 '+ 


A\x) = 




( 2 . 2 ) 


(2.3) 


We define ( V e,+ , W e ' + ) for x < -L as the unique solution of 12.21 with initial values at x = —L given by 
the values, (V e + , W e ' + ) , of the singular solution 11.12) at x = —L. Furthermore, we define U B ’ + for x < —L 
by the second equation in 12.1) . U e ’ + = ^-W e,+ . In this way the singular solution U 6, + is defined for all 

itl. 


Assuming that uj > c\9\ the system 12. II has two independent plane-wave solutions for x < —L 

,2 




= i,±* 


U! 


c 2 k(oj) 


s ±lfc( “ )x , k{u>):=\l^ T -9 2 . 


(2.4) 


Note that the condition w > c\9\ implies fc(w) G M which yields truly propagative plane-waves. Defining 
U± = 737 - W± . we obtain a basis of plane-wave solutions to 12 . II for x < —L, 


Z ± = (u±,V±',Ul). (2.5) 

Then, for some functions / (x), R(co) 

U 9 ’ + = I(u) Z e + + R(u) Z® , for x < -L. (2.6) 

The functions I(uj),R(uj) can be easily computed in terms of the values at x = —L of of U e + in II. 121 . we 
omit the details. 

To consider a wave packet with frequency support in an interval [wo, ui\ we make the following assump¬ 
tions to make sure that the upper-hybrid resonance takes place at a unique point x^ for all ui G [wo,wi]. 
We have that 


a := 


UJ 2 (UJ + L0 h ) 

c 2 (w 2 — w 2 ) 


(w 


Wft). 


Recall that ou c := oo p := and that toh := ^Jlo 2 + w 2 is the upper-hybrid frequency. As above, 

assume that the density of electrons N e is zero for x < — L, that it is strictly increasing for —L < x < H , that 
it is constant for x > H and that it is twice continuously differentiable. We take w 0 > !jJ c and wi < u h (H). 
Then, for w G [w 0 ,wi], a u {-L) > 0 and (//) < 0 with a unique point x w G (- L,H) where a u (x u ) = 0 and 


a'(x u ) := ,) < 0 . 
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We apply the results in Section 1 and 1151 for each w e [tu 0 , cui] by translating x u to zero and we make 
explicit the dependence in ui of the singular solution denoting it by U® ,+ . We assume that w 0 > c\8\ and we 
consider a wave packet that is incident from the left. 


P v (t) := [ e~ iut A(u) (U«>+, tp) du, (2.7) 

Jr 

where A(u) is an integrable function with support contained in [wo,wi], (this means that A(u) can be 
different from zero only in [wo,wi] )• The function tp is continuously differentiable on R, with values in 
C 3 , and has compact (closed and bounded) support. By (-, •) we denote the duality between the space of 
distributions and the test functions, 


(U e '+,^)= / U e j+(x)tp(x)dx. 

J R 

Note that by 12.61 A(u>) I (ui) is the amplitude of the plane wave that is incoming from —oo and that A(lo) R(lo) 
is the amplitude of the reflected plane wave. Furthermore, as U e + is a distribution it does not make sense 
to evaluate it in point-wise sense nor to compute its L 2 norm. This is the reason why we introduced in 
(22) the test function tp. As tp has support in a compact set, say K, in intuitive terms by testing the wave 
packet with tp we obtain information in how the wave packet behaves in AT as a function of time. Note 
that, in a naive way we could try to define, 

[ e~ iut A(oj)\J e j + (x)dx, (2.8) 

Jr 

but, of course, since U 0 + is a distribution this integral has to be given a meaning in mathematical sense. 
The quantity P v (t) in 12.71 can be understood as a proper definition (or as as a regularized version) of 
12.81 . where by changing the function tp we test the wave packet in different compact set in R. Actually 
as Uf/+ is a function away from the hybrid resonance at x = 0, we could study the behavior of the wave 
packet in point-wise sense, or in local L 2 norms, away from the upper-hybrid resonance, but this is not 
the point, because we wish to consider precisely the influence of the upper-hybrid resonance, and hence 
we do not dwell on this issue here. 

We will use the Riemann-Lebesgue theorem. For the reader’s convenience we recall it here. Let tp be a 
continuously differentiable function defined on R, with values on C 3 and with compact support. Let ip be 
its Fourier transform, 

ip(t) := f e~ ltx ip(x) dx. 

Jr 

Integrating by parts one proves that for some constant C, 

\m\<C^, feR. 
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Furthermore, if ij) is integrable, approximating it in the L 1 norm by continuously differentiable functions 
with compact support one proves that, 

lim = 0. (2.9) 

£->-±oo 

Let us decompose Uf/ + between a regular part and part uf /2 that contains all the singularity 

U®’ + = U*;j + U®’+ (2.10) 

where, denoting by xo( x ) the characteristic function of the set O, i. e. , Xo( x ) = 1 , x G O, xo( x ) = 0, x ^ O, 

U®’,t “ U^’ + (*)X(- 0 o,-L ) (®) + K’ + , v e j+, w 0 j + ) X [-l,oo )(*), (2.11) 


and 


U e ’+ — 

9 • — 


:= ( p - k mF + rfE7 fo(a: " ) - 0 ' 0 J 


We have proven in Il5l that. 


l( 


0 ,+ y0,+ 

a UJ 5 U UJ 5 W La 


) \\»<c, 


( 2 . 12 ) 


(2.13) 


where the constant C is uniform for u; in compact sets. This is the reason why is referred to as the 
regular part. Then, 

P v (t) = P vA (t) + P v , 2 (t) (2-14) 

where, 

P v ,i(t) := f e- iut A{u) (v e j+, v )du, (2.15) 

P v At)~ f e- iut A{u)(u e j+,<p)duj. (2.16) 

Then, using also 12.61 we obtain that is bounded for w G [wo,uu] and since A (A) is integrable, it 

follows from 12.91 that. 


lim P v , i(f) = 0. 

£—>■±00 


(2.17) 


We now consider P v , 2 {t) where we have to take into account the singularity of the solution. 

For e > 0 let p w {±e) satisfy: a(p u (±e)) = ±e. Note that if e is small enough, Voj G [wo,wi], p u (±e) is 
unique and there is a p 0 > 0 such that, \p u {±e)\ > p 0 . Moreover, a(x) is one to one for x G \pu(e), Pu>{— e)] 


and 


ol ' ( x ) 


is uniformly bounded for w G [wo,wi], x G \p u {e), p u (—e)]. Observe that this holds because N e is 


continuously differentiable and p u (e) < Pu,(x u ) < p u ( s) and a'(x u ) < 0 . 
We decompose P v p{t) in the following way 

P^2(t) = PS(t)+Pg(t), 


(2.18) 
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with. 


P${t) := 

with r] > e, and 


-rr (x(-l,p^W)(^) + X(pU-v ),°°)( x )) <p(x) + ™ ip{x u ) ) dx 

a[x) v ' ' OL r {x u ) ' 


P, 


?l(t) := [ e lujt A(w) h(u)dui, where h{uS) := lim 
’ J eto 


rPm{e) y 


' PuW) 


a(x) 


f 

tp(x)dx + j 
J o, 


e iut A(uj)duj, 

Pu(-v) i 


Pu(-e) 


a(x) 


ip(x)dx 


(2.19) 


( 2 . 20 ) 


Since the function inside the brackets in the right-hand side of 12.191 is uniformly bounded for u € 
[wo,wi] and A(ui) is integrable, it follows from 12.91 that. 


lim P£ 2 (f) = °- 


( 2 . 21 ) 


Let fiu{z) be the inverse function to a u (x) i.e., P 0J {a U] (x)) = x, x G \ftu{y), p u (—r))\. Then, changing the 
variable of integration into 2 := a u (x) we obtain that the function h(/x) in 12.201 is equal to 


h(u>) = — lim 

ej.0 


f g 1 v(PUz)) 
!-n za'(Pu(z)) 


dz 


: 1 f <p(Pu(z)) _ 


n 1 tp(/3 u (z)) 
e za'{P u {z)) 

1 


dz 


efo [J z \a’(P u (z)) a'(J) u (0)) 


dz 


r* 1 {<p(j3 u (z)) 


and then. 


r 1 / 1 1 


z \a'03u,{z)) a'(f3u,(0)) 

^CM 0 ))) dz. 


<p(/3o;(0))) dz 


'(&,(*)) a '(p u ( 0 )) 

Hence, h(u) is bounded for u G [w 0 , uq] and as A(u) is integrable, by 12.91 . 

I™ P^lit)=0. 

t—>± OO 

Then, we have proven, 

THEOREM 2.1. For every continuously differentiable ip with compact support. 


( 2 . 22 ) 


(2.23) 


lim P v (t) = 0. 

t->±oo ^ 


(2.24) 


Proof: The theorem follows from 12.141 . 12.171 . 12.181 . 12.211 and 12.231 . 

Equation 12.241 . with appropriate functions tp, implies that that as I —> ±00 the part of the wave packet 
in [-R, -L] tends to zero for every R> 0. Since there is no heating for x < —L this means that the reflected 
part of the wave packet actually travels to —00 in the x direction. Furthermore, the wave packet tends to 
zero, as t —> ±00 for x > L. Since we assume that our problem is coercive (non propagative) for x >11 
the wave packet cannot travel to plus infinite in the x direction. This implies that for large times the 
transmitted part of the wave packet is dissipated by collisions. 

Note that the condition w 0 > c\6\ is only used in order that for x < —L we have an incoming and a 
reflected plane wave, see 12.61 . Equation 12.241 also holds for w 0 < c|0| with a similar proof. We omit the 
details. 
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